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Electronic surface states in one-dimensional two-band TBA model are studied by use of the Green 
function method. The local density of states (LDOS) at successive atoms in a semi-infinite chain, 
even in the case of atoms distant from the surface, is found to be clearly different from that observed 
in an unperturbed (infinite) chaini. The surface atom occupancy is calculated self-consistently^, 
with the effect of electron-electron interactions taken into account. The electron-electron interactions 
are shown to have a significant impact on the conditions of surface state existence. 



I. INTRODUCTION 

The presence of electronic surface states has a substan- 
tial effect on the properties of solids. Many apparently 
surprising features of mesoscopic and nanoscopic systems 
- in which surface effects are particularly conspicuous, the 
surface representing a significant part of the whole - can 
be elucidated by the conditions of existence of surface 
states. 

Adsorption and reconstruction processes render the de- 
scription of real surfaces difficult. The properties of real 
systems can often be precisely reproduced through nu- 
merical simulations, which, however, do not provide the 
explanation of the mechanism of generation of surface 
states and their effect on bulk states. 

One of the basic models used for description of elec- 
tronic properties of solids is the tight binding approx- 
imation (TBA) model. The pioneering studies on the 
conditions of electronic surface state existence, based on 
a single-band model of finite crystal, were reported by 
Goodwin^. The single-band model was then generalized 
by Artmani, who introduced a double-band model to 
investigate the existence of two types of surface states: 
Shockley states, which are induced only by breaking 
the translational symmetry of the crystal, and Tamm 
states, generated as a result of introducing an additional 
perturbation^. A breakthrough was marked by the pa- 
per by Kalkstein and Soven^, in which Green's function 
formalism was used for the determination of properties 
of surface and bulk states in a semi-infinite crystal with 
perturbed surface. 

The problem of multi-electron effects and their impact 
on surface state existence conditions has not yet been 
exhaustively discussed. Papers on chemisorption, which 
is a related issue, are available, though^ii. As in the 
case of chemisorption, the simples way of including the 
impact of the multi-electron effects on the conditions of 
surface state existence consists in incorporating interac- 
tion of electron with charge density into the Coulomb 
model. Introducing the Hartree potential into the TBA 
model of a semi-infinite crystal amounts (in the simplest 
case) to surface atom site energy renormalization. 

This study is focused on surface states in a ID semi- 
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FIG. 1: The model of semi-infinite ID crystal with two-atom 
unit cell. The surface site (n = 0) is occupied by an adatom 
with site energy q'^. The alternating resonance integral sign 
is a consequence of interaction between orbitals s and p. 



infinite atom chain being a model of ionic crystal with 
two atoms in the unit cellS.. Multi-electron effects are 
taken into account in Hartree approximation only. The 
surface atom occupancy and site energy values are found 
through self-consistent calculations using Green's func- 
tion formalism^. 



II. MODEL 

The model assumes non-zero resonance integral values 
for neighboring sites only. Orbitals s and p are alter- 
nately attributed to successive chain sites. Hence, reso- 
nance integrals for successive site pairs alternate in sign, 
taking values /? or —(3. Site energy associated with s or 
p orbital is denoted ae or cto, respectively (cf. FigOJ. 

A wave function in the TBA model is assumed to be a 
linear combination of atomic functions: 

N 

|fc)=^(c2„|2n)-|-C2„+i|2n-H)), (1) 

n 

the sum involving all the two- atom unit cells. Expressed 
in the atomic function basis, the Hamiltonian of an infi- 
nite (unperturbed) chain has the following form: 

Ha = ^(ae|2n)(2n|+ao|2n+l)(2n + l| + 

n 

-K/3|2n)(2n+l|-/3|2n+l)(2n|). (2) 
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A surface introduced into the system is regarded as a 
perturbation breaking the infinite chain into two separate 
semi-infinite ones^: 

H = Ho + V (3) 
V ^ (a^-ae)(|-l)(-l| + |0)(0|) + 

+/?(|-1)(0| + |0)(-1|). (4) 

The non-zero values of elements V (0, 0) and V (—1, —1) 
allow for adsorption of atom of different type. 

Derived from the secular equation for Hamiltonian Hq , 
the expansion coefficients c„ and the dispersion relation 
read as follows: 



C2rn = Ae'^'l^ 



X = ±^Jt'^ + 2 - 2cos( 



(5) 
(6) 



where 



A=J^^ i3 = ^l^^, (7) 



X + T 



and 



X, \X\>T 

r, |X|<T 



(8) 



Parameters X and 9 represent dimensionless energy and 
wave vector, respectively: 
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a and r are defined as follows: 



ae + ao Q!e — Q-o 

a = : , T 



2 ' ■ 2/3 
The Greenian matrix elements: 



E-E{k) 



(9) 



(10) 



(11) 



expressed in the atomic function basis, read within the 
energy bands: 



Go(2TO,2n) 



X + T t^~" 



Go(2m+ l,2n+ 1) = ' < 

Go(2m,2n+ 1) 



13 t^-ty' 

,n- 



f3 t<-U 



(12) 
(13) 
(14) 



Go(2m + l,2n) = ~ < < ^ , (15) 



where 



i> = Z ± sign(X)iyZ2~T, 



(16) 



and 



Z = 



COSf 



(17) 



Selected Greenian matrix elements for perturbed 
(semi-infinite) crystal can be found from Dyson's equa- 
tion: 



G = Go + GqVG. 
The diagonal elements read: 



(18) 



G(m, m) = Go(to, m) + (19) 
Go(0, m) (Go(m, 0)^(0, 0) + Go(m, -1)F(-1, 0)) 
1 - Go(0, 0)F(0, 0) - Go(0, -1)F(-1, 0) 

This allows the determination of the local density of 
states (LDOS): 



p{X, m) = — 7r9 [/3G(m, to)] 



(20) 



and the surface state occupancy in successive chain sites: 
lahcucha: 



{n{m)) — Res [G(to, to), Xg] , 



(21) 



where Xg is the surface state energy determined from the 
condition of G(to, m) zeroing. 

Self-consistent renormaHzation of site energy at succes- 
sive sites is necessary for electron-electron interactions to 
be taken into account. By defining 



f3 



"5 ^rn — *-^ei ^o; Q^e; Q^oi ^e; ■ 



we get: 



= Trn + U {n{m)} 



(22) 



(23) 



where J7 is a parameter defining interaction of electron 
with charge density. Surface state localization is equiva- 
lent to state occupancy fading inward the crystal. There- 
fore, the highest site energy gradient is expected at the 
surface. In the first approximation, site energy modifica- 
tion can concern the surface atom only. 

The surface perturbation parameter can be expressed 
as follows: 

Ae ^ "^-"^ =To-T2„, m= 1,2,3,.... (24) 

When interaction of electron with charge density is taken 
into account: 



K^Tl^-T2m, m= 1,2,3,. 



(25) 



In the case considered here (r = 1), the surface state 
energy is expressed by the following formula: 



1 + Af ± yi - 4A^ -I- 6Af + iA'jTAf 
Xs± = , (26) 

Xs+ and Xs- being the solutions vaHd for 
|A^ + 1/2| > y/E/2 and A^ < 0, respectively. 
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FIG. 2: Surface state levels around the upper band (gray 
area). The solid and dotted lines represent the state levels 
found with or without the electron-electron interactions taken 
into account, respectively. The dashed lines delimit the region 
in which no surface states exist. 



III. RESULTS 

Computations were performed at r = 1 {ae — Oo = 2/3). 
FiglSl shows the surface state energy, Xg, plotted versus 
the surface perturbation. No Shockley states are found 
to exist in the model discussed'^, as no surface states are 
found in the absence of perturbation. Tamm states, in- 
duced through modifying the surface atom site energy, 
Ag, are found to emerge from the upper energy band 
(associated with ae ). 

The surface states emerging from the bottom edge of 
the band are induced by arbitrarily small perturbation 
value. For surface states to be induced above the upper 
band, however, the perturbation value must be positive 
and fulfill the condition Ag > 1/2(^/5— 1). Perturbation 
values from the interval < A'g < l/2(-\/5— 1) correspond 
to surface state non-existence. The solid and dotted lines 
in Fig in represent surface state levels found with or with- 
out the electron-electron interactions taken into account, 
respectively. Clearly, the multi-electron effects (in the 
Hartree approximation) boost the surface state energy 
levels, resulting in weakened or strengthened localization 
of surface states below or above the upper band, respec- 
tively. However, the electron-electron interactions have 
no effect on the interval of surface perturbation parame- 
ter values at which surface states are found to exist. 



The effect of the surface on the electronic states in the 
considered chain is the most evident in the LDOS spec- 
trum. FigEI shows the LDOS plots obtained for the four 
sites closest to the surface {n = ... 3) . Three different 
perturbation values are assumed, corresponding to sur- 
face state appearing below the upper band (Ae — —0.75), 
not induced at all (Ag = 0.5), and induced above the 
upper band (Ag — 1). The solid and dotted lines rep- 
resent the LDOS calculated with multi-electron effects 
taken into account or neglected, respectively. As a re- 
sult of including the multi-electron effects, the LDOS in 
the upper band is increased; at the same time, the oc- 
cupancy of the surface state below the band decreases, 
and the state moves towards the band edge. An oppo- 
site effect is found to occur for the surface state above 
the upper band: the LDOS in the band is found to de- 
crease, while the surface state occupancy increases and 
the state moves inwards the gap. Comparing relative oc- 
cupancy changes at successive sites, one notes stronger 
localization in states closer to the band edge. Because of 
band asymmetry, the surface state occupancy should be 
compared between sites either even or odd. With multi- 
electron effects taken into account, the occupancy ratio 
of site (n = 0) (the surface atom) to site n = 2 is found 
to increase or decrease for states above or below the up- 
per band, respectively. Note that even in the absence of 
surface states the presence of the surface still affects the 
LDOS spectrum. The van Hove singularities at the band 
edges are eliminated, and LDOS minima appear inside 
the bands. Even for sites very distant from the surface, 
the LDOS differs significantly from that in an infinite 
(unperturbed) chain^. 



IV. CONCLUSIONS 

Focused on the effect of electron-electron interac- 
tions on the conditions of surface state existence (in 
the Hartree approximation), the above-presented study 
shows that the localization of surface states generated 
above or below the upper energy band can be increased 
or decreased, respectively, by the electron-electron inter- 
actions in the considered model. However, the interval of 
surface perturbation parameter values corresponding to 
surface state existence is found to remain unaffected by 
taking these interactions into account in the calculations. 
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FIG. 3: The local density of states, (7rp(X)), at successive chain sites starting from the surface (n = 0, 1, 2, 3). Results obtained 
at three different perturbation values, Ae, are grouped into columns. The dashed line represents the LDOS for infinite chain. 
The dotted and solid lines represent the LDOS found for semi-infinite chain with electron-electron interactions neglected or 
taken into account, respectively. The arrows indicate the surface state occupancy. 



